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The N/Z dependence of projectile fragmentation at rclativistic energies has been studied 
in a recent experiment at the GSI laboratory with the ALADiN forward spectrometer 
coupled to the LAND neutron detector. Besides a primary beam of ^^^Sn, also secondary 
beams of ^■^^La and ^'^'^Sn delivered by the FRS fragment separator have been used in 
order to extend the range of isotopic compositions of the produced spectator sources. 
With the achieved mass resolution of AA/A x 1.5%, lighter isotopes with atomic num- 
bers Z < 10 are individually resolved. The presently ongoing analyses of the measured 
isotope yields focus on isoscaling and its relation to the properties of hot fragments at 
freeze-out and on the derivation of chemical freeze-out temperatures which are found 
to be independent of the isotopic composition of the studied systems. The latter result 
is at variance with the predictions for limiting temperatures as obtained with finite- 
temperature Hartree-Fock calculations. 

* E-mail: w.trautmann@gsi.de 
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1. Introduction 



The present interest in isotopic effects in nuclear multifragmentation is partly moti- 
vated by the importance of the nuclear equation of state for astrophysical processes. 
Supernova simulations or neutron star models require inputs for the nuclear equa- 
tion of state at extreme values of density and asymmetry for which the predictions 
differ widely.^"^ Fragmentation reactions are of interest here because they per- 
mit the production of nuclear systems with subnuclear densities and temperatures 
which, e.g., largely overlap with those expected for the explosion stages of core- 
collapse supernovae.'^ Laboratory studies of the properties of nuclear matter in the 
hot environment, similar to the astrophysical situation, are thus becoming feasible, 
and an active search for suitable observables is presently underway. Of particular 
interest is the density-dependent strength of the symmetry term which is essen- 
tial for the description of neutron-rich objects up to the extremes encountered in 
neutron stars (refs.'^'^~^ and contributions to this workshop). 

Here, new results from a recent experiment performed with the ALADIN spec- 
trometer at the GSI laboratory will be discussed in which the possibility of using 
secondary beams for reaction studies at relativistic energies has been explored. 
Beams of ^"^Sn, ^^^Sn, ^^"'La, and -'^^^Au as well as Sn and Au targets were used 
to investigate the mass and isospin dependence of projectile fragmentation at 600 
MeV per nucleon. The neutron-poor radioactive projectiles ^°^Sn and ^^''La were 
produced at the Fragment Separator FRS by fragmentation of a primary beam of 
^^^Nd and delivered to the ALADIN experiment. 

The results presented in the following, after a brief discussion of global observ- 
ables, will be restricted to quantities deduced from the measured yields of resolved 
isotopes. They will, in particular, include isoscaling and its relation to the proper- 
ties of hot fragments at the low-density freeze-out, as well as chemical temperatures 
deduced from double ratios of isotopic yields and the consequences for the limiting- 
temperature concept in multifragmentation. Preliminary results from the present 
experiments have been presented previously. ^"^^^^ 

2. N/Z dependence of global fragmentation observables 

The dependence on the isotopic composition of the system is rather small for global 
observables of the studied reactions. This is shown in Fig.[T]for the multiplicity of the 
produced intermediate-mass fragments and for Zmax, both as a function of Abound- 
Here 2^max denotes the largest atomic number Z within a partition while the sorting 
variable .Abound = SZ^ with Z^ > 2 is related to the impact parameter and inversely 
correlated with the degree of excitation of the produced spectator system. The 
multiplicities exhibit the universal rise and fall of fragment production - (Fig. [H 
left panel), and only a slightly steeper slope in the rise section distinguishes the 
neutron-rich case of -'^^'^Sn from the other two systems. As confirmed by statistical 
model calcul atio ns, this difference is related to the evaporation properties of excited 
heavy nuclei.'^'^ Neutron emission prevails for neutron-rich nuclei which leads to a 
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Fig. 1. Mean multiplicity <Mimf> of intermediate-mass fragments 3 < Z < 20 produced in the 
fragmentation of ^"^'^■^^Sn and ^^*La at 600 A MeV as a function of Abound (l^ft panel). The 
right panel shows the correlations of the mean Z of the largest fragment with .Z^ound with both 
quantities being normalized with respect to the projectile Z. 

concentration of the residue channels, with smaU associated fragment multiphcities, 
in a somewhat narrower range of large 2^bound than in the case of the neutron-poor 
systems. The evaporation of protons reduces .Abound since protons are not counted 
therein. 

This effect is, nevertheless, small and nearly invisible in the correlation of 
<-^max> with Abound (Fig. [TJ right panel). There, the transition from predomi- 
nantly residue production to multifragmentation becomes apparent as a reduction 
of <Zniax> with respect to Abound which occurs between ^bound/^^proj = 0.6 and 
0.8. The range of isotopic compositions of the studied nuclei N/Z = 1.14 to 1.48 is 
not large enough for significant variations of the reaction mechanism to appear. A 
shrinking of the coexistence zone in the temperature-density plane is predicted for 
neutron-rich matter but observable consequences can only be ex pected for asym- 
metries far beyond those presently available for laboratory studies! ^ ^ This fact, 
on the other hand, has the important consequence that the basic reaction process is 
the same for all the present systems, a prerequisite for the interpretation of isoscal- 
ing and its relation with the symmetry energy. This point will be further quantified 
in the section on temperatures. 



3. N/Z dependence of light fragment production 

The mass resolution obtained for projectile fragments entering into the acceptance 
of the ALADIN spectrometer is about 3% for fragments with Z < 3 and decreases 
to 1.5% for Z > 6. Masses are thus individually resolved for fragments with atomic 
number Z < 10. The elements are resolved over the full range of atomic numbers 
up to the projectile Z with a resolution of AZ < 0.2 obtained with the TP-MUSIC 
IV detector 
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Fig. 2. Inclusive mean values <N>/ Z of light fragments with 3 < Z <9 produced in the fragmen- 
tation of io^'i24sn (600 A MeV, ALADIN), SBpe and ^ssy (both 1 A GeV, FRS, from refs.lMSl) 
as a function of the fragment Z (left panel). The right panel shows the results for Z = 8 as a 
function of the N/Z value of the projectile. The lines represent the trend of the data (full line) 
and <N>/Z = (N/Z)p^aj (dashed). 



The inclusive mean neutron-to-proton ratio <N>/Z of light fragments in the 
range 3 < ^ < 9 is presented in Fig. [51 The values for Z — 4 have been corrected 
for the missing yield of unstable ^Be by including an estimate for it obtained from 
a smooth interpolation over the identified yields of ^^^~^^Be. This correction has a 
negligible effect for the case of ^"^Sn because <N>/Z of the detected beryllium iso- 
topes is close to 1 already. For ^■^"'Sn the value of <N>/ Z for Z = 4 is lowered from 
1.24 to 1.18 which makes the systematic odd-even variation more clearly visible in 
the neutron rich case. The odd-even staggering is, however, much more pronounced 
in the case of the neutron-poor ^"^^Sn. The strongly bound even-even nuclei attract 
a large fraction of the product yields during the secondary evaporation stage.L^^ 
This effect is, apparently, larger if the hot fragments are already close to symmetry, 
as it is expected for the fragmentation of ^°^Sn.LLiJ 

Inclusive dat a " ^ ° l obtained with the FRS fragment separator at GSI for 238 u 
and ^^Fe fragmentations on titanium targets at 1 A GeV bombarding energy confirm 
that the observed patterns are very systematic (Fig. [21 left panel). Nuclear structure 
effects characteristic for the isotopes produced combine with significant memory 
effects of the isotopic composition of the excited system by which they are emitted. 
This has the consequence that, because of its strong variation with Z, averaged 
neutron-to-proton ratios <N>/ Z for unspecified ranges of Z are not very useful 
for studying nuclear matter properties. For this purpose, techniques will have to be 
used which cause the nuclear structure effects to cancel out. Selecting a particular 
element, e.g. Z = 8, is already sufficient to reveal the linear correlation of <N>/Z 
with the N/Z of the projectile (Fig. [21 right panel). According to the Statistical 
Fragmentation Model (SMIVP^, the correlation should be much stronger for the 
excited fragments at the breakup stage and, in fact, not deviate much from the 
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dashed <N>/Z = (7V/Z)proj line in the figureM^The reduction is due to sequential 
decay and its tendency of directing the isotope distributions closer to the valley of 
stability. For a quantitative analysis, a precise modeling of these secondary processes 
will, therefore, be required. 
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Fig. 3. Scaled isotopic ratios S{N) for Li to O isotopes from the fragmentation of -^^^Sn and ^"^Sn 
projectiles at 600 A MeV (ALADIN, left panel) and for H to B isotopes from l^Q -|- Il2,i24gj^ ^^ 
E/A = 300 and 600 MeV (INDRA, right panel, from refl^^^) as a function of the neutron number 
N . The centrality selection was obtained with .Abound (ALADIN, left panel, five bins as indicated) 
and according to charged particle multiplicity as measured with the full detector (INDRA, right 
panel, three bins with "central" indicating 6/6max < 0.4). The dashed lines in both panels are the 
results of exponential fits according to Eq. JTJ. Only statistical errors are displayed; note also the 
offset factors of multiples of three in the right panel. 



4. Isoscaling 

The phenomenon of isoscaling has been shown to be common to many different 
types of heavy-ion reactions. ^'^^^'^ It is observed by comparing product yields Yi 
from reactions which differ only in the isotopic composition of the projectiles or 
targets or both. Isoscaling refers to an exponential dependence of the measured 
yield ratios R2i{N,Z) on the neutron number iV and proton number Z of the 
detected products. The scaling expression 



R2i{N, Z) = Y2{N, Z)/Yi{N, Z) = C ■ exp{aN + [3Z) 



(1) 



describes rather well th e m easured ratios over a wide range of complex parti- 
cles and light fragments.'^ The accuracy with which the isoscaling relation is 
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obeyed is conveniently judged by regarding the scaled isotopic ratios S{N) = 
R2iiN, Z)/exp{pZ), obtained by dividing out the common Z dependence after 
fitting according to eq. H]). Isoscaling means that all the scaled yield ratios will 
follow the same exponential dependence on the neutron number N. The results for 
the pair of reactions with ^^"^Sn and ^°'^Sn projectiles, sorted into five bins of Abound 
are shown in Fig. [3l left panel. Alternatively, one may also regard the behavior of 
the exponential slopes a{Z) of the yield ratios for individual elements or /3(iV) for 
individual chains of isotones. As demonstrated in Fig. |4] for the inclusive yields 
from the same pair of reactions, the deviations from the common mean value are 
impressively small. Besides the precision with which isoscaling is observed in the 
range of fragments from Li to Ne, also the monotonic decrease of the slope of S{N) 
with centrality represents a remarkable observation. Qualitatively, it resembles that 
observed for the reactions ^^C + ^^^'^^"^Sn at 300 and 600 MeV per nucleon, studied 
with INDRA at GSI {rei^, right panels of Fig.[3l). 
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Fig. 4. Isoscaling parameters a{Z) and f3{N) obtained by individually fitting the A'^ and Z depen- 
dences of isotopes and isotones, respectively, for inclusive yields from the fragmentation of ^^''Sn 
and ^'^''Sn projectiles at 600 A MeV. The dashed lines represent the mean values for isotopes with 
Z from 3 to 10 and isotones with A'^ from 3 to 12. Only statistical errors are displayed. 



The statistical model offers a simple physical explanation for the appearance 
of isoscaling in finite systems. Charge distributions of fragments with fixed mass 
numbers A, as well as mass distributions for fixed Z, are approximately Gaussian 
with average values and variances which are connected with the temperature, the 
symmetry-term coefficient, and other parameters. (Here the symmetry-energy term, 
Esym — liA — 2Z)'^/A with the coefficient 7, is that of the liquid-drop description 
of the nascent fragments at freeze-out and 7 = 25 MeV is normally used.^i^ The 
mean values of the fragment distributions depend on the total mass and charge of 
the systems, e.g. via the chemical potentials in the grand-canonical approximation, 
while the variances depend mainly on the physical conditions reached, the temper- 
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ature, the density and possibly other variables. For example, the charge variance 
cz ~ vI^^/^t) obtained for fragments with a given mass number A is only a func- 
tion of the temperature and of the symmetry-term coefficient since the Coulomb 

contribution is very small!^ This relation of isoscaling with the symmetry energy 

l91 99 9"^ 97 9^ 
has attracted considerable interest recently. ^^ ^^ ^' ^o 

In the grand-canonical approximation, assuming that the temperature T is 
about the same (see below), the scaling parameters a and (3 are given by the dif- 
ferences of the neutron and proton chemical potentials for the two systems divided 
by the temperature, a — Afi^/T and /3 — A/ip/T. The proportionality of A/i,i 
and thus of the isoscaling parameters with the coefficient 7 of the symmetry-energy 

term has been obtained from the statistical int erp retation of isoscaling within the 

91 i94l 

SMM and Expanding-Emitting-Source Model^-^and confirmed by an analysis of 

reaction dynamics.l^ The relation is 

AMn = Mn,2 - A^„,i « 47(|| - ||) = 47A(ZVA2) (2) 

where Zi,Ai and .^2,^2 are the charges and mass numbers of the neutron-poor 
and neutron-rich systems, respectively, at breakup. The difference of the chemical 
potentials depends essentially only on the coefficient 7 of the symmetry term and 
on the isotopic compositions. A proportionality of a and 1/T, with T derived from 
double isotope ratios, has been observed for light-ion (p, d, a) induced reactions 
at bombarding energies in the GeV range.l^ Also the symmetry-term coefficient 
7 Ri 22.5 MeV, deduced according to eq. ([2]), supported the statistical approach. 
These data were inclusive and a value close to the standard value was to be expected 
because the fragment multiplicities are small in this case.^^ 

The present reactions studied with ALADIN cover the full rise and fall of frag- 
ment production while the ^^C induced reactions studied with INDRA cover mainly 

the rise up to .^bound_~ ■^proj/2 with multiple fragment production being dominant 

11 01 
in central collisions.'''^ In both cases, the measured temperatures rise slightly with 

increasing centrality but not sufficiently fast in order to compensate for the decrease 

of the isoscaling parameter a. The product A/in — a-T, therefore, decreases. With 

the assumption that the isotopic compositions remain close to those of the original 

ii2,i24gj-^ targets, the evolution of the apparent symmetry-term coefficient 7app, i.e. 

without corrections for the effects of sequential decay, was followed into t he regime 

of multifragmentation with the ^^C induced data measured with INDRA.I^It was 

again found to be close to the normal-density coefficient 7 « 25 MeV for peripheral 

collisions, associated with small fragment multiplicities, but dropped to lower values 

of 15 to 17 MeV for the central impact parameters leading to multifragmentation. 

This result is confirmed by the preliminary analysis of the new data for projectile 

fragmentation measured with ALADIN. The product a ■ T drops by about 20% 

through the rise of fragment production but then remains fairly constant once the 

regime of multifragmentation has been reached (^bound/^proj < 0.5). 
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Fig. 5. Apparent temperatures TucLi (left panel) and TboLI (right panel) as a function of Zbound 
for the three reaction systems produced with ^"^'^^^Sn and ^^^La projectiles. Only statistical 
errors are displayed. 



5. N/Z dependence of the nuclear caloric curve 

Two examples of double-isotope temperatures deduced from the measured isotope 
yields are shown in Fig. [5] as a function of .Abound- Besides the frequently used 
^HeLi (left panel), determined from ■^■"'He and ^'^Li yields also the results for TbcLI 
are displayed (right panel). For TbcLi the isotope pairs of ^^^Be and ^^^Li are used 
which each differ by two neutrons. The double difference of their binding energies 
amounts to 11.3 MeV and is nearly as large as the 13.3 MeV in the TneLi case. 
The apparent temperatures are displayed, i.e. no corrections for secondary decays 
feeding the ground states of these nuclei are applied. 

Both temperature observables show the same smooth rise with increasing cen- 
trality that was observed earlier in a study of ^^^Au fragmentations.'^ According 
to TbolIj the temperatures at the chemical freeze-out are identical for all three 
reaction systems and have the same dependence on .Abound- This is not equally 
visible in TneLi which exhibits slightly larger values for the neutron rich case of 
^^''Sn. Inspection of the single isotope ratios shows, however, that the ^'''He ratio 
represents an exception in that it varies much less than the other ratios with the 
neutron-richness of the system. A possible explanation is offered by the observa- 
tions shown in Fig. [21 The strong population of even-A'^-even-Z nuclei, evident from 
the behavior of <N>/Z, suggests that also "^He nuclei are abundantly produced. 
The effect is stronger in the neutron-poor cases for which the apparent TuoLi will 
be accordingly lower. Taking this nuclear-structure effect into account is likely to 
reduce the differences between systems also for this thermometer. 

The observed invariance of the breakup temperature with the isotopic composi- 
tion of the system is of interest because it is opposite to what is expected according 
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to the finite-temperature Hartree-Fock calculations of Besprosvany and Levitl^ 
The limiting temperatures obtained for the stability of excited nuclei are strongly 
dependent on the Coulomb pressure generated by the protons they contain. Along 
isobars, the limiting temperatures decrease rapidly toward the proton rich side and, 
along the valley of stability, they decrease with increasing mass because the effect of 
the increasing Z dominates over that of the decreasing Z/A. For the nuclei studied 
here, a difference of about 2 MeV is predicted between the limiting temperature of 
^^"'Sn and those for the proton rich ^°^Sn and ^^''La or for the heavier ^^^Au!^ A 
difference of this magnitude is clearly not seen in the data (Fig. [5|) . Neither does 
a comparison with the apparent temperatures measured for ^^^Au fragmentations 
(Fig. 7 in ref."^^ provide evidence for the predicted mass dependence. The tem- 
peratures are the same, within errors, when regarded as a function of the scaled 
quantity .^bound/^proj ■ This, however, is suggested by the Abound scaling of fragmen- 
tation observables if a com parison of sources of different mass but equal excitation 
is intended to be made.L!^ 



6. Discussion 

It is one of the basic assumptions of isotopic reaction studies that the mechanism as 
such remains unchanged and that effects related to the asymmetry dependence of 
the nuclear forces can be isolated by changing nothing but the isotopic composition 
of the system. This view is supported by the observed invariance of the chemical 
breakup temperatures, here for the case of spectator fragmentation. The interpre- 
tation of the breakup temperatures as a manifestation of the limiting temperatures 
predicted by the Hartree-Fock model,'^on the other hand, is not equally supported. 
This interpretation was derived from the 2'bound dependence of the temperature 
which sugg ested a mass dependence because of the changing mass of the spectator 



lOOl 

system.l^^ Later on, the predicted mass dependence was found to be confirmed in 
a compilation of chemical break-up temperatures from several experiments,'^' in- 
cluding the earlier ALADIN data, which led to their use as a general reference for 
testing predictions of nuclear stability with respect to temperature (see ref.'^- for 
references). The new data for the A w 120 region and their comparison with ^^^Au 
fragmentations suggest that the observed variation with decreasing .Zbound is more 
related to the increasing excitation energy than to the decreasing mass of the pro- 
duced spectator system. Breakup temperatures fairly independent of the isotopic 
composition or mass of the system are, e.g., predicted by the statistical fragmenta- 
tion models which are based on the concept of a phase-space driven instability.'^ 
The discussion of the isoscaling results has focussed on mainly three topics, 
among them primarily the effects of sequential decay for the isoscaling parameters 
and for the symmetry term, as well as the role of the surface term in the symmetry 
energy. Surface terms gain in importance as the mean mass of the produced frag- 
ments decreases in multi-fragment decays. Thirdly, also a possible modification of 
the source N/Z due to preequilibrium emissions prior to breakup will affect and 
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most probably reduce the isoscaling parameters. So far, however, none of these 
mechanisms has been identified as providing an obvious explanation for the large 
drop that is o bserved. Calculations with the microcanonical Markov-chain version of 
the SMMjI^even show that sequential decay will cause the isoscaling coefficients to 
increase if the symmetry term governing the mass distributions of the hot fragments 
is reduced.'^ Transport models predict that the isotopic composition at breakup 
should not deviate by more than a few percent from the original value.l^ New 
experimental results, furthermore, indicate that the isotopic dependence of the sur- 
face term is reduced, with respect to that needed to describe ground-state ma sses , 
as the excitation of the system increases.!^ This is also expected theoretically.'^ 

A decrease of the isoscaling parameter a with increasing violence of the colli- 
sion, beyond that expected from the simultaneous increase of the te mperature, has 
also been observed for reactions at intermediate energy (refs! ' ' and references 
therein) and a systematics is emerging. A very recent interpretation of these re- 
sults arrives at the conclusion that the reduced density, rather than the eleva ted 
temperature, causes the symmetry term to be smaller than the standard value.'^ 
This analysis uses a temperature independent potential part of the symmetry term 
and a kinetic part in the form of a Fermi gas, and it is argued that the isoscal- 
ing phenomenon itself is evidence for the fact that properties of the considered 
homogeneous nuclear matter are observable in frag mentation reactions. 

In the Statistical Multifragmentation Model, L!^ normal-density fragments statis- 
tically distributed within an expanded volume are considered. Here a reduction of 
the symmetry term in their liquid-drop description can be imagined to be caused by 
fragment modifications in the hot environment, including deformations or the effects 
of nuclear interactions between them. The neglect of such effects in the actual codes 
should probably be regarded as an idealization as the assumption of equilibrium at 
the final freeze-out stage requires nuclear interactions in order to achieve it. The 
consequences of a reduced symmetry term for predicted global fragment observ- 
ables, as e.g. multiplicity, are rather small for this class of models. The partitioning 
is predominantly driven by the surface term in the fragment descripti on while a 
variation of the symmetry term affects mainly the isotopic distributions! ' l ^'^ l '^" ! 

Further studies will be needed in order to put some of the assumptions on a 
firmer basis and to quantitatively establish the reduction of the symmetry term 
in multi-fragment breakups. The sequential decay corrections are obviously impor- 
tant but also the evolution of the isotopic compositions during the reaction process 
deserves attention. It will have to be explored whether an experimental reconstruc- 
tion of the neutron-to- proton ratio of the detected spectator systems will be feasible 
using the coincident neutron data measured with the LAND detector.!^ 

This work has been supported by the European Community under contract 
No. HPm-CT-1999-00001, by the Polish State Committee for Scientific Research 
(KBN) under contract No. 2P03B11023, and by the Polish Ministry of Science and 
Higher Education under Contracts No. 1 P03B 105 28 (2005 - 2006) and N202 160 
32/4308 (2007-2009). 
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